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Summary: 

We  determined  the  equation  of  state  for  ballistic  gelatin  using  the  Brillouin 
scattering  spectroscopy  with  a  diamond  anvil  cell  by  measuring  the  pressure  dependence 
of  the  sound  velocities  from  ambient  pressure  to  12  GPa  at  room  temperature.  We 
extended  these  measurements  to  a  temperature  range  of  0  to  100  °C  between  ambient  and 
12  GPa.  We  analyzed  the  Brillouin  data  using  a  high  temperature  Vinet  equation  of  state 
and  obtained  the  bulk  modulus,  its  pressure  derivative,  thermal  expansion  coefficient  and 
the  temperature  derivative  of  bulk  modulus  at  ambient  pressure.  We  also  applied  the  Mie- 
Gruneisen  method  to  the  P,  V,  T  data  and  obtained  Gruneisen  parameter,  Debye 
temperature  and  its  thermal  pressures  for  the  gelatin. 
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I.  Introduction 


The  ubiquity  of  polymer-like  materials  in  science  and  society  is  a  clear  testament 
to  the  importance  of  these  unique  materials.  The  diversity  of  polymeric-material 
applications,  which  range  from  optics  and  electronics  to  garden  hoses  and  soda  bottles, 
provides  a  practical  impetus  for  fundamental  research  in  polymer  physics  and  chemistry. 
While  the  chemistry  of  polymeric  chain  formation  is  primarily  dictated  by  the  reactivity 
of  the  monomer  units,  the  resultant  physical  properties  of  a  given  polymer  are  highly 
reliant  on  network  topology  and  presence  or  absence  of  filler  materials.  In  addition, 
external  stimuli,  such  as  variable  temperatures  and  pressure,  have  pronounced  influence 
on  their  mechanical  behavior.  Examples  include  thermal  aging  phenomena,  pressure- 
induced  crystallization,  shock-compressed  polymorphism,  and  thermoplasticity.  Most  of 
the  biomaterials  (such  as  gelatin)  can  be  analogous  to  polymers. 

The  term  “gelatin”  is  used  to  refer  to  a  yellow-white  powder,  100%  protein, 
extracted  by  acid  or  base  hydrolysis  from  collagen,  the  main  matrix  material  of  animal 
skin,  bone,  and  connective  tissue.  The  term  also  refers  to  a  water  solution  of  this  powder. 
Gelatin  is  an  aggregate,  meaning  there  is  no  single  molecular  size  and  formula.  Instead, 
gelatin  is  a  collection  of  large  protein  molecules  similar  to  amino  acid  composition.  In  a 
given  gelatin  sample,  the  molecular  weight  ranges  between  17,000  and  300,000  Daltons. 
Thus  only  the  average  molecular  weight  can  be  used  in  characterization  of  the  gelatin. 

Gelatin,  one  of  the  most  versatile  hydrocolloids,  has  a  wide  variety  of  uses  in 
many  industries.  It  is  used  in  the  food  industry  for  making  edibles  such  as  desserts, 
marshmallows,  jellies,  and  aspics.  It  is  used  by  the  pharmaceutical  industry  to 
manufacture  capsules,  tablets,  and  bandages.  Gelatin  is  also  an  important  raw  material  in 
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the  photographic  industry.  Another  use  for  gelatin  is  as  the  “model”  of  human  tissue  in 
bullet  penetration  experiments  (Ref.  1).  This  is  the  so-called  “ballistic  gelatin.”  It  is 
necessary  to  obtain  the  equation  of  state  (EOS)  for  ballistic  gelatin  to  describe  its  elastic 
behavior  under  pressure. 

Since  sound  waves  are  intimately  associated  with  changes  in  pressure  and  density, 
Brillouin  scattering  measurements  can  provide  precise  pressure- volume  data  for  a  variety 
of  materials.  The  work  presented  here  supplies  a  necessary  extension  to  the  currently 
limited  EOS  data  available  for  gelatin.  In-situ  sound  velocities  are  directly  related  to 
elastic  constants  which,  in  turn,  arise  directly  from  the  curvature  of  an  intermolecular 
potential.  Given  these  relationships,  sound  velocity  measurements  can  provide  useful 
insight  into  the  microscopic  structure  of  a  material.  The  EOS  and  mechanical  properties 
of  gelatin  -  have  been  determined  from  ambient  pressure  to  approximately  12  GPa  using 
a  diamond-anvil  cell  (DAC)  in  conjunction  with  Brillouin  scattering.  Besides  elastic 
properties  and  P-V  EOS  data,  Brillouin  scattering  can  also  provide  insight  into  the 
dynamic  relaxation  processes  in  gelatin  (Ref.  2). 

II.  Experimental  details 

(1)  Diamond  anvil  cell 

Figure  1  shows  the  simple  schematic  diagram  of  DAC  used  in  this  experiment.  It 
consists  of  three  major  components:  the  diamond  anvils,  a  piston,  and  cylinder  parts.  The 
diamond  anvils  are  separately  fixed  to  the  piston  and  cylinder  parts.  A  metal  gasket  with 
a  hole  in  its  center  is  used  as  a  sample  chamber.  When  the  piston  and  cylinder  are  pressed 
toward  each  other,  the  gasket  as  well  as  the  pressure  medium  and  samples  will  be 
compressed,  which  will  generate  pressure  on  the  sample.  In  this  experiment,  we  did  not 
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use  a  pressure  medium  because  gelatin  itself  is  very  soft  and  would  play  the  role  of 
pressure  medium.  We  used  a  standard  ruby  fluorescence  method  to  determine  the 
pressure.  The  samples  of  gelatin  (20%)  are  prepared  with  the  methods  described  in  Ref. 
1.  Gelatin  samples  are  loaded  into  a  DAC  (Fig.  1)  with  no  pressure  medium  added.  Pre¬ 
indented  (12  GPa)  stainless  gaskets  with  gasket  hole  of  about  150  microns  are  used. 
Figure  1(b)  shows  some  details  of  diamond  anvils,  gasket,  and  thermocouples. 

(2)  Brillouin  scattering  system 

Brillouin  light  scattering  is  generally  referred  to  as  inelastic  scattering  of  an 
incident  optical  wave  field  by  thermally  excited  elastic  waves  in  a  sample.  Light 
scattering  from  thermally  excited  acoustic  waves  was  theoretically  predicted  by  L. 
Brillouin  (Ref.  3).  Later,  Gross  (Ref.  4)  experimentally  confirmed  the  prediction  in 
liquids  and  in  crystals.  Brillouin  spectroscopy  did  not  progress  until  the  1960s,  due  to  a 
lack  of  intense  light  sources.  In  1960,  the  invention  of  the  laser  solved  the  light  source 
problem;  since  then,  Brillouin  light  scattering  techniques  have  earned  reasonable 
attention  and  have  been  widely  used  in  many  research  fields. 

In  the  case  of  a  transparent  solid  (crystals  or  glasses),  the  incident  laser  light  will 
induce  dynamic  fluctuations  in  the  strain  field  to  bring  about  fluctuations  in  the  dielectric 
constant,  and  these  in  turn  translate  into  fluctuations  in  the  refractive  index  due  to  the 
elasto-optic  scattering  mechanism.  These  fluctuating  optical  inhomogeneities  result  in 
inelastic  scattering  of  the  light  as  it  passes  through  the  solid.  The  phonons  present  inside  a 
solid  move  in  thermal  equilibrium  with  very  small  amplitudes,  creating  fluctuations  in  the 
dielectric  constant,  which  is  viewed  as  a  moving  diffraction  grating  by  an  incident  light 
wave.  Therefore  Brillouin  scattering  can  be  explained  as  a  combination  of  Bragg 
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reflection  and  Doppler  effect.  The  moving  grating  scatters  the  incident  light  giving  rise  to 
scattered  photons  with  Doppler  shifted  frequencies  Av.  Brillouin  spectrum  gives 
frequency  shift  (Av)  of  the  thermal  phonon,  and  its  wavelength  ( d  space)  can  be 
determined  from  the  experiment  geometry.  For  instance,  with  n;  as  the  refractive  index  in 
the  direction  of  incident  light,  ns  as  the  refractive  index  in  the  scattering  direction,  and  n 
as  the  velocity  of  acoustic  phonons,  the  Brillouin  shift  Av  is  given  as  follows: 

=  cos60  C1) 

A 

where  2  is  the  wavelength  of  incident  light  and  6  is  the  scattering  angle. 

The  simple  schematic  arrangement  for  a  Micro-Brillouin  scattering  system  is 
shown  in  Figure  2,  and  some  experimental  details  can  be  found  in  Refs. 5- 10.  This  is  a 
non-contact  measurement  technique  that  exploits  light  scattering  to  probe  the  elastic 
properties  of  liquids  and  solids  (crystals).  Light  from  a  single  mode  Ar-ion  laser  (2=514.5 
nm)  was  used  as  the  excitation  source  with  the  average  power  less  than  100  mW.  The 
sample  was  placed  symmetrically  with  respect  to  the  incoming  and  collected  light  such 
that  the  difference  vector  of  the  incoming  and  detected  light  was  in  the  plane  of  the 
sample.  The  scattered  light  was  collected  by  lens  L2,  analyzed  by  a  3+3  tandem  Fabry- 
Perot  interferometer,  detected  by  a  photon  counting  photomultiplier,  and  was  outputted  to 
a  multichannel  scalar  (Fig.  2).  Spectra  were  taken  at  each  pressure  and  a  temperature 
point,  and  collected  for  10  to  120  minutes. 

Figure  3  exhibits  the  symmetric  scattering  geometry  employed  in  these 
experiments,  the  scattering  angle  was  80°  degrees.  Since  samples  were  placed 
symmetrically  with  respect  to  the  incoming  and  collected  light,  and  for  this  particular 
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geometry  the  frequency  shift  of  the  incident  light  is  independent  of  the  refractive  index  of 
the  samples,  the  Eq.  (1)  can  be  reduced  to 

Av  =  (2u//l)sin(#/2),  (2) 

where  n  is  sound  velocity,  2=514.5  nm  (laser  wavelength),  and  (9=80°  is  the  scattering 
angle. 

The  temperature  up  to  100  °C  was  generated  by  a  ribbon  heater  which  wrapped 
around  DAC  as  shown  in  Fig.  1(a).  The  temperature  was  measured  by  a  thermocouple 
which  was  placed  between  the  diamond  anvil  and  the  gasket  (close  to  the  sample)  as 
shown  in  Fig.  1(b).  Since  the  heater  was  uniformly  wrapped  over  the  DAC  body,  the 
temperature  could  be  considered  as  uniformly  distributed.  Since  diamond  is  an  excellent 
heat  conductor,  the  measured  temperature  can  be  considered  close  enough  to  the  sample 
temperature.  The  fluctuations  of  temperature  were  below  ±  3  °C  at  100  °C  during  the 
measurement. 

III.  Results 

The  microscopic  photo  images  taken  at  pressure  0,  2.1,  and  11  GPa,  respectively 
are  shown  in  Fig.  4.  The  images  at  0  and  1 1  GPa  show  that  the  gelatin  clearly  is  isotropic. 
However,  the  image  at  2.1  GPa  shows  that  the  gelatin  starts  to  crystallize  at  this  pressure. 
The  images  at  2.1  and  11  GPa  were  taken  during  different  runs  of  experiments.  The 
resulting  images  can  be  understood  as  follows:  the  samples  are  compressed  to  around  2 
GPa  and  stay  at  this  pressure  regime,  the  samples  will  crystallize;  if  they  are  compressed 
very  quickly  and  pass  through  this  pressure  regime,  the  samples  will  not  crystallize  and 
will  go  to  a  glassy  state. 
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A  typical  set  of  Brillouin  spectra  is  shown  in  Figure  5.  There  are  two  pairs  of 
peaks  in  each  spectrum  and  an  additional  Rayleigh  peak  at  zero  frequency  (elastic 
scattering).  One  pair  of  peaks  corresponds  to  the  longitudinal  acoustic  mode  (L-mode), 
whereas  the  other  pair  corresponds  to  the  transverse  acoustic  mode  (T-mode). 

In  additions,  figure  6  exhibits  the  pressure  dependence  of  Brillouin  frequency 
shifts  at  room  temperature.  Figure  7  shows  the  pressure  dependence  of  sound  velocities  at 
selected  temperatures  between  273  and  373  K.  Figure  8  shows  the  temperature 
dependence  of  sound  velocities  at  selected  pressures.  Now,  we  have  the  pressures, 
temperatures,  and  sound  velocity  data  (calculated  from  Eq.  2)  between  ambient  and  12 
GPa,  273  and  373  K.  The  P-V-T  EOS  fitting  will  be  discussed  in  detail  in  the  next 
section. 

IV  Discussion 


1 1]  is  the  ratio  of  specific  heat  at  constant  pressure  and  constant  volume. 

To  calculate  the  density  changes  from  0  to  100  °C  between  0  and  12  GPa,  we 
have  to  solve  the  following  two  problems  first:  (1)  the  initial  density  for  each  temperature 
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between  0  and  100  °C  at  ambient  pressure;  (2)  the  initial  velocities  for  each  temperature 
between  0  and  100  °C  at  ambient  pressure. 

We  solved  these  problems  as  follow:  (1)  the  initial  density  was  measured  at  30  °C 
at  ambient  pressure,  p=0.991  +/-  0.006  g/cm3.  Our  initial  density  is  very  close  to  the  data 
from  Ref.  12.  We  applied  the  T-dependence  of  density  from  Ref.  12  and  extrapolated  it 
up  to  100  °C  at  ambient  pressure.  The  result  is  shown  in  Fig.  9  with  the  inset  table.  (2) 
The  initial  sound  velocity  is  1340  m/s  at  30  °C  at  ambient  pressure.  Since  gelatin  is  not 
stable  above  40  °C  (Ref.l),  it  is  difficult  to  measure  the  sound  velocity  above  40  °C  at 
ambient  pressure.  Instead  we  measured  the  temperature  dependence  of  sound  velocity 
from  0  to  100  °C  at  0.53  GPa  as  shown  in  figure  10.  We  assumed  that  the  T-dependence 
of  sound  velocity  at  ambient  pressure  should  be  similar  to  the  one  obtained  at  0.53  GPa 
due  to  the  small  pressure  difference.  Therefore,  by  using  the  T-dependence  of  sound 
velocity  at  0.53  GPa  and  combined  with  the  initial  sound  velocity  measured  at  30  °C  at 
ambient  pressure,  we  extrapolated  the  T-dependence  of  sound  velocity  from  0  to  100  °C 
at  ambient  pressure.  The  results  are  shown  in  Fig.  10  (b).  Consequently,  we  calculated  the 
pressure  dependence  of  density  (1/V)  at  each  selected  temperature.  The  results  are  shown 
in  Figure  11  and  table  I.  There  are  very  small  difference  in  density  between  -2  and  100 
°C.  In  the  next  two  sections,  we  will  discuss  how  we  analyze  the  P-V-T  Brillouin  data 
using  high  temperature  Vinet  EOS  and  the  Mie-Gruneisen  method. 

(2)  P-V-T  EOS 

The  P-V-T  EOS  are  needed  in  order  to  discuss  the  intermolecular  interaction  and 
potential  for  simple  fluid  or  solid  materials.  However,  in  the  case  of  gelatin,  the  EOS  can 
only  give  us  the  average  information  about  intermolecular  interaction  and  potential.  Here 
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we  analyze  the  P-V-T  data  with  two  methods:  one  extends  the  room-T  Yinet  EOS  to  a 
higher  temperature  by  considering  the  thermal  expansion,  the  other  is  based  on  the  idea 
of  thermal  pressure,  which  is  that  the  total  pressure  on  the  sample  is  the  sum  of  static 
pressure  and  pressure  due  to  thermal  expansions. 

(a)  High  temperature  Vinet  EOS 

Yinet  equation  of  state  is  expressed  as  follow: 

\~(VIVVn  1 

P  =  3K0[  ;  °2'3  ]exp[-«-l)(l-(F/F0)1/3)]  (4) 

where  Ko  is  the  bulk  modulus,  Ko  ’  is  its  derivative  at  zero  pressure,  and  Vo  is  the  initial 
volume  at  room  temperature.  We  fitted  the  pressure  dependence  of  density  at  room 
temperature  by  using  the  Vinet  EOS  with  three  parameters  and  obtained  the  bulk 
modulus  of  4.9  GPa,  its  derivative  at  zero  pressure  of  8.2,  and  the  initial  volume  of  1.004 
(cm3/g).  For  the  temperature  effect,  the  high  temperature  Vinet  EOS  incorporates 
temperature  effect  into  the  third  order  Vinet  EOS,  and  is  expressed  as  the  following: 

]-(v/v  y/3  t. 

P  =  3K,rl  (v,v  V  ]exp[-«r  -1)(1  (5) 

where  Vot,  Kot,  and  Fot  are  isothermal  bulk  modulus,  its  pressure  derivatives,  and  initial 
volume  at  temperature  T  and  ambient  pressure,  respectively.  Assuming  that  the  second 
and  third  order  pressure  derivatives  of  the  bulk  modulus  are  negligible,  then  Kot  and  Kot 
are  given  by  the  properties  at  ambient  conditions 


K„=K„+(^)p(T-302) 

oT 

(6) 

Kqt  =  ^0 

(V) 

11 


where  Kq  and  Ko  are  the  isothermal  values  at  ambient  conditions.  The  temperature 
derivatives  of  the  bulk  modulus  are  assumed  to  be  constant  throughout  the  whole 
temperature  range.  Zero-pressure  volume  Vot  is  expressed  as 

Vor  =  V0  exp  j^aTc/r  (8) 

(Zj,  =  oCq  +  cc^T  +  0C2T  2  (9) 

where  To  is  the  initial  volume  at  ambient  conditions,  and  aj  is  the  thermal  expansion 
coefficient  at  ambient  pressure,  and  expressed  by  constant  parameters  ao,  a\,  and  a.2. 
Since  the  experiment  was  carried  out  in  a  relatively  narrow  temperature  range  from  0  to 
100  °C,  the  coefficients  a\  and  «2  can  be  set  to  zero.  The  combination  of  Eqs.  (5)  to  (9)  is 
fitted  to  the  P-V-T data.  A  least-squares  fit,  using  program  EOS-fit  Version5.2  (Ref.  13), 
yielded:  K0= 4.6  +/-  0.2  GPa,  Kd ’=8.7  +/-  0.2,  T0=  1.004  +/-  0.002  cm3/g,  a0=(40.4  +4.0) 

xlO-5  K'1,  and  (— )P  =-0.009  +/-  0.002  GPaK'1.  We  separately  fitted  the  Vinet  EOS  to 
dT 

pressure  dependence  of  density  at  each  temperature.  We  obtained  (Tot,  Ko t,  Aot’)  for 
each  temperature;  the  results  were  summarized  in  Table  II  &  III. 

(b)  Mie-Gruneisen  model 

According  to  the  Mie-Gruneisen  theory,  the  total  pressure  Ptot(V,  T)  at  a  given 
volume  and  temperature  can  be  expressed  as  a  sum  of  static  compression  at  a  reference 
temperature  Ps t(T,  To)  and  the  thermal  pressure  APth(T,  T)  at  a  given  volume  V: 

Ptot  (V,  T)  =  Pst  (V,  T0  )  +  A Pth  0 V,T )  (10) 

For  the  static  pressure  term,  Pst,  we  use  the  Vinet  fit. 

The  Debye  model  can  be  used  to  describe  thermal  pressure  term  in  (10).  Indeed,  it 
has  been  suggested  that  the  Debye  model  is  a  good  approximation  for  solids  such  as  MgO 
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or  MgSiC>3  or  many  other  minerals.  We  only  can  assume  that  the  Debye  model  may  also 
be  a  possible  approximation  for  gelatin  due  to  its  elastically  isotropic  properties.  With 
this  model,  the  first  terms  in  (10)  is: 

1-(V/V)l/3 

P  =  3 K0 [  ;  °2; 3  ] exp[—  (K0  - 1)(1  ~(V /F0)1/3)]  (11) 

V  '  V o)  ^ 

and  the  second  term  is 


A Pth(V,T) 


riV) 

v 


[Eth(V,T)-Eth(V,T0)] 


(12) 


where  Kq  is  the  isothermal  bulk  modulus,  Kq  is  its  pressure  derivative,  and  y  is  the 
Gruneisen  parameter.  The  subscript  zero  denotes  ambient  conditions.  The  vibrational 
energy  can  be  calculated  from  the  Debye  model 


9 nRT  ten 

~  ( JJrf  J0  e"  -\  * 


(13) 


The  Debye  temperature  is  given  by  6,  n  is  the  number  of  atoms  in  a  formula  unit,  and  R  is 
the  gas  constant.  The  Gruneisen  parameter  y  is  assumed  to  be  a  function  of  volume  only 

r  =  ro(^r  (14) 

where  q  is  assumed  to  be  constant.  The  Debye  temperature  can  also  be  described  as  a 
function  of  volume 


e=eatxf(-  —)  (is) 

q 

In  this  approach,  the  fitting  parameters  are:  Vo,  Ko,  Kq\  y0,  q,  and  6*o-  This  EOS  has  the 
advantage  that  it  yields  higher  order  thermoelastic  terms  in  an  internally  consistent  way. 
The  first  three  terms  ( Vo ,  Ko,  Ko ’)  are  obtained  from  300  K  data.  Although  the  EOS  is  not 
sensitive  to  the  choice  of  6o,  calculating  the  Debye  temperature  for  gelatin  is  an  important 
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step.  In  general,  Debye  temperature  is  calculated  for  crystalline  structures  (metal, 
crystals,  and  minerals  [Ref.  11])  using  the  lattice  vibrational  properties.  However,  the 
gelatin  does  not  have  crystalline  structure.  We  may  be  able  to  estimate  the  value  for 
gelatin  as  follows:  Since  we  measured  the  sound  velocity  at  ambient  pressure  and  room 
temperature  and  the  longitudinal  velocity  ul=1340  m/s,  shear  velocity  Dt  =0  m/s.  Dt  can 
be  replaced  by  dl  and  using  the  following  equation 


v 


3  t,3  + 


v, 


(16) 


we  obtained  the  Debye  velocity  Dd=1933  m/s.  The  Debye  temperature  can  be  expressed 
as  follows: 


0  K  At:  M/p  D 


(17) 


where  k,  h,  and  N  are  Boltzman's  constant,  Plank's  constant,  M  is  the  molecular  mass,  p  is 
the  number  of  atoms  in  the  molecular  formula,  and  Avagadro's  number,  respectively. 
Since  gelatin  does  not  have  crystalline  structure,  we  use  the  initial  volume  at  ambient 
conditions.  The  average  weight  of  gelatin  molecular  is  175000  Dalton  which  equals  to 
3*10'  kg.  The  number  of  molecules  in  initial  volume  is  3.33*10  .  Using  those 
numbers,  we  estimated  the  Debye  temperature  as:  100  K<  Go  <  200  K. 

The  combination  of  Eqs  (7)-(12)  were  fitted  to  the  P-V-T data.  Three  different  0o 
(100,  150,  and  200  K)  were  chosen.  The  fitting  results  remained  unchanged. 
Consequently,  Go  is  set  to  150  K.  This  allows  us  to  obtain:  Ko=A.9  +/-  0.2  GPa,  Kq,= 8.2 
+/-  0.2,  Fo=1.004  +/-  0.002  cm3/g,  Gruneisen  parameter  y=0.1 7,  and  q=l  (fixed).  The 
calculated  pressures  as  well  as  the  thermal  pressures  are  shown  in  Table  Y.  Figure  12 
shows  the  thermal  pressure  as  a  function  of  temperature.  For  example,  the  thermal 
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pressure  at  372  K  and  9  GPa  is  about  0.6  GPa.  This  number  is  within  the  ±3c  of 
measured  error  in  pressure  (a=0.2  GPa).  This  indicates  that  the  temperature  range  in  the 
current  study  is  relatively  narrow,  in  order  to  well  constrain  the  Gruneisen  parameter  and 
its  q,  it  is  necessary  to  extend  the  temperature  range  up  to  or  above  600  °C  . 

(c)  Elastic  properties 

The  isotropic  material  has  two  independent  elastic  constants:  Cn  and  Cn-  In 
general: 

pvl=Cn  (18) 

PU2r=\(.C„-Cu)  (19) 

where  p  is  density,  ul  is  the  longitudinal  sound  velocity,  and  ut  is  the  transverse  sound 
velocity.  The  Poisson’s  ratio  can  be  expressed  as  follows: 

o_l  -  2 vj  (20) 

2  (o2l-o2t) 

We  calculated  Cn,  C12,  and  o  as  the  function  of  pressure  at  room  temperature  shown  in 
Fig.  13.  They  both  increase  linearly  with  pressure;  Cn  increases  with  pressure  faster  than 
C12.  On  the  other  hand,  the  Poisson’s  ratio  decreases  with  pressure.  The  value  of 
Poisson’s  ratio  is  0.36  at  10  GPa.  It  is  comparable  with  that  of  silicate  glass  at  the  same 
pressure  (Ref.  5). 

Conclusion 

We  determined  the  EOS  for  ballistic  gelatin  using  the  Brillouin  scattering 
spectroscopy  using  a  diamond  anvil  cell  to  measure  the  pressure  dependence  of  the  sound 
velocities  from  ambient  pressure  to  12  GPa  (in  a  80°  symmetric  scattering  geometry)  at 
room  temperature.  We  calculated  the  pressure  dependence  of  elastic  constants  and 
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Poisson’s  ratio  at  room  temperature.  We  extend  this  method  to  a  temperature  range  of  0 
to  100  °C  between  ambient  pressure  and  12  GPa.  We  analyzed  the  Brillouin  data  by 
using  a  high  temperature  Yinet  EOS  fitting  to  P-V-T  data  and  obtained  the  thermal 
expansion  coefficient  (first-order  approximation).  We  also  used  the  Mie-Gruneisen 
method  and  obtained  the  Gruneisen  parameter  y  for  gelatin. 
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(a) 


Scattering  light 


Figure  1.  (a)  Simple  diagram  of  diamond  anvile  cell.  The  Brillouin  scattering  geometry  is 
80  degrees:  symmetric  scattering,  (b)  Diamond  anvils. 
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Figure  2.  Simple  diagram  of  Brillouin  scattering  system.  Ar+  ion  laser  is  operated  at 
514.5  nm  and  power  is  less  than  100  mW. 
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Figure  3.  A  diagram  of  the  symmetric  scattering  geometry.  The  incident  light  enters  the 
sample  with  the  angle  612  ,  and  by  Snell’s  law  the  following  equation  will  hold: 
n o s i n( /9/2)=n s i n( r/j/2 ) ,  where  no  is  the  refractive  index  of  air  and  n  is  the  refractive  index  of 
samples.  Therefore  the  Brillouin  shifts  will  have  the  following  form: 
2/7  &  2u 

= - sin(^>/2)  =  — sin (91 2)  =  Arr  The  Brillouin  shifts  can  be  defined  uniquely 

2  2 

by  the  scattering  angle  between  incident  and  scattering  light. 
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150  um 


Figure  4.  Photo  image  of  gelatin  in  a  DAC  cell  at  selected  pressure. 
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Intensity  (arb.  units) 


Figure  5.  A  typical  Brillouin  spectrum.  The  elastic  scattering  line  at  the  zero  frequency. 
The  Brillouin  lines  at  around  5  and  12.5  GHz  for  transverse  and  longitudinal  acoustic 
mode,  respectively. 
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Frequency  shift  (GHz) 


Figure  6.  Pressure  dependence  of  the  Brillouin  frequency  shifts.  Solid  and  open  circles 
correspond  to  the  longitudinal  and  transverse  acoustic  mode,  respectively. 
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Temperature  (°C) 


Figure  7.  The  temperature  dependence  of  Brillouin  shifts  which  correspond  to  the 
longitudinal  mode  at  selected  pressure.  Symbols  are  experiment  data  and  lines  are  the 
linear  fits. 


Figure  8.  Pressure  dependence  of  sound  velocity  at  selected  temperatures.  Symbols  are 
the  experiment  data  and  lines  are  the  guide  to  eyes. 
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20  40  60  80  100 

Temperature  (°C) 


Figure  9.  Showing  the  T-dependence  of  density  at  ambient  pressure.  The  data  (open 
circle)  are  from  Ref.  12.  The  dashed  line  is  the  extrapolation  from  0  to  100  °C. 
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Figure  10.  (a)  The  typical  Brillouin  spectra  as  the  function  of  temperature  at  0.53  GPa; 
(b)  it  shows  the  T-dependence  of  Brillouin  frequency  shifts.  Symbols  are  data  and  solid 
line  is  the  linear  fit;  (c)  T-dependence  of  linewidth  at  0.53  GPa.  The  behavior  of 
linewidth  relates  to  the  relaxation  process  in  the  sample  (Ref.  2).  It  is  beyond  the  current 
study. 
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Figure  13.  Pressure  dependences  of  elastic  constants  and  Poisson’s  ratio  at  room 
temperature. 


Table  I.  Showing  pressure,  temperature,  calculated  density,  volume,  and  bulk  sound 
velocity,  L-mode  and  T-mode  velocities.  The  pressure  has  an  averaged  error  of  ±  0.2 
GPa,  the  Brillouin  has  the  error  of  50  m/s  in  measurements. 


Pressure 

Temperature 

Density 

Volume  (cm3 

Bulk  velocity 

L-Mode 

T-Mode 

(GPa) 

(K) 

(g/cm3) 

/g) 

(m/s) 

(m/s) 

(m/s) 

0 

271 

1.008 

0.992 

1434 

1434 

0 

0.05 

271 

1.026 

0.975 

2029 

2029 

0 

0.53 

271 

1.103 

0.907 

3562.45 

3562.45 

0 

1.42 

271 

1.169 

0.855 

3818.84 

4479.19 

2027.23 

3.23 

271 

1.275 

0.784 

4563.02 

5314.3 

2359.12 

5.47 

271 

1.370 

0.73 

5214.68 

6041.32 

2641.68 

6.53 

271 

1.407 

0.711 

5518.34 

6387.11 

2785.2 

7.59 

271 

1.440 

0.694 

5758.35 

6714.97 

2991.51 

8.8 

271 

1.476 

0.678 

5905.13 

6872.84 

3045.33 

10.48 

271 

1.522 

0.657 

6128.82 

7137.9 

3168.66 

0 

302 

0.995 

1.005 

1363 

1363 

0 

0.05 

302 

1.015 

0.985 

1953 

1953 

0 

0.53 

302 

1.098 

0.911 

3446.29 

3446.29 

0 

1.42 

302 

1.165 

0.858 

3835.69 

4325.35 

1731.22 

3.23 

302 

1.270 

0.787 

4561.45 

5153.28 

2076.56 

5.47 

302 

1.366 

0.732 

5169.93 

5957.45 

2563.64 

6.53 

302 

1.404 

0.712 

5456.99 

6274.54 

2682.04 

7.59 

302 

1.438 

0.695 

5695.81 

6624.54 

2929.16 

0° 

oo 

302 

1.473 

0.679 

6062.53 

6964.78 

2969.08 

10.48 

302 

1.518 

0.659 

6199.21 

7288.15 

3318.91 

0 

320 

0.989 

1.011 

1322 

1322 

0 

31 


6.53 

372 

1390 

0.719 

5369.97 

6236.86 

2747.07 

7.59 

372 

1424 

0.  702 

5800.99 

6778.21 

3036.36 

8.8 

372 

1458 

0.  686 

6193.81 

7145.53 

3085.69 

10.48 

372 

1505 

0.  664 

6126.47 

7127.14 

3153.86 

Table  II:  Bulk  modulus,  its  derivative,  initial  volume,  thermal  expansion  coefficient,  and 
the  temperature  derivative  of  bulk  modulus  at  ambient  pressure  are  obtained  via  high- 
temperature  Yinet  EOS  fit. 


271  K 

302  K 

320  K 

333  K 

372  K 

Yot  (GPa) 

5.42  ±0.4 

5.2  (fix) 

5.9  ±0.6 

8.1  ±0.7 

3.4  ±0.4 

Vot’ 

8.1  ±0.4 

8.3  (fix) 

7.8  ±0.4 

6.8  ±0.4 

9.3  ±0.6 

E0T(cm3/g) 

0.991  ±0.003 

0.996(fix) 

0.991  ±0.008 

0.965  ±0.005 

1.039  ±0.005 

a0(10'5  K'1) 

32.3  ±4.6 

(Xi 

0 

(dKT/dT)P 

-0.008±0.002 

Table  III:  The  initial  density  is  set  to  fixed  number,  the  Vinet  EOS  fit  yields  the  bulk 
modulus  and  its  derivative. 


271  K 

302  K 

320  K 

333  K 

372  K 

K0j  (GPa) 

5.4  ±0.3 

4.6  ±0.3 

4.7  ±0.3 

3.5  ±0.4 

3.1  ±0.3 

Yot’ 

8.1  ±0.3 

8.7  ±0.5 

8.4  ±0.4 

9.8  ±0.6 

9.6  ±0.5 

E0T(cm3/g) 

0.992  (fix) 

1.004  ±0.002 

1.011  (fix) 

1.015  (fix) 

1.046  (fix) 

a0(10'5  K'1) 

40.4  ±  4.6 

ai 

0 

(dKr/dT)P 

-0.009±0.002 

Table  IV:  High  temperature  Vinet  EOS  fitted  to  the  P-V-T  data.  The  recalculated 
pressures  were  shown  in  table.  Fitting  parameter:  Ko=4.6  +/-  0.2  GPa,  V)’ =8. 7  +/-  0.2, 

F0=  1.004  +/-  0.002  cm3/g,  a0=(40.4  +4.0)  xKf5  K1,  and  {^-)P  =-0.009  +/-  0.002 

dT 

GPaK'1. 

Reference  temperature  =  302.00  K,  Pressure  in  GPa. 


V  (cm3/g)  T  (K)  P  (obs)  P  (calc)  AP 


0.9920 

271.0 

0.0000 

-0.0013 

0.0013 

0.9747 

271.0 

0.0500 

0.0913 

-0.0413 

0.9066 

271.0 

0.5300 

0.6429 

-0.1129 
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0.8554 

271.0 

1.4200 

1.3522 

0.0678 

0.7843 

271.0 

3.2300 

3.0666 

0.1634 

0.7299 

271.0 

5.4700 

5.3473 

0.1227 

0.7107 

271.0 

6.5300 

6.4595 

0.0705 

0.6944 

271.0 

7.5900 

7.5687 

0.0213 

0.6775 

271.0 

8.8000 

8.9129 

-0.1129 

0.6570 

271.0 

10.4800 

10.8473 

-0.3673 

1.0050 

302.0 

0.0000 

-0.0036 

0.0036 

0.9852 

302.0 

0.0500 

0.0952 

-0.0452 

0.9107 

302.0 

0.5300 

0.6810 

-0.1510 

0.8584 

302.0 

1.4200 

1.4021 

0.0179 

0.7874 

302.0 

3.2300 

3.0961 

0.1339 

0.7321 

302.0 

5.4700 

5.3889 

0.0811 

0.7123 

302.0 

6.5300 

6.5271 

0.0029 

0.6954 

302.0 

7.5900 

7.6686 

-0.0786 

0.6789 

302.0 

8.8000 

8.9723 

-0.1723 

0.6588 

302.0 

10.4800 

10.8541 

-0.3741 

1.0111 

320.0 

0.0000 

0.0019 

-0.0019 

0.9901 

320.0 

0.0500 

0.1040 

-0.0540 

0.9124 

320.0 

0.5300 

0.7088 

-0.1788 

0.8584 

320.0 

1.4200 

1.4565 

-0.0365 

0.7880 

320.0 

3.2300 

3.1435 

0.0865 

0.7358 

320.0 

5.4700 

5.2703 

0.1997 

0.7163 

320.0 

6.5300 

6.3551 

0.1749 

0.6993 

320.0 

7.5900 

7.4706 

0.1194 

0.6826 

320.0 

8.8000 

8.7460 

0.0540 

0.6623 

320.0 

10.4800 

10.5860 

-0.1060 

1.0152 

332.0 

0.0000 

0.0053 

-0.0053 

0.9930 

332.0 

0.0500 

0.1110 

-0.0610 

0.9124 

332.0 

0.5300 

0.7382 

-0.2082 

0.8576 

332.0 

1.4200 

1.5040 

-0.0840 

0.7868 

332.0 

3.2300 

3.2251 

0.0049 

0.7326 

332.0 

5.4700 

5.4834 

-0.0134 

0.7128 

332.0 

6.5300 

6.6217 

-0.0917 

0.6954 

332.0 

7.5900 

7.7974 

-0.2074 

0.6794 

332.0 

8.8000 

9.0632 

-0.2632 

0.6553 

332.0 

10.4800 

11.3422 

-0.8622 

1.0460 

372.0 

0.0000 

-0.0466 

0.0466 

1.0194 

372.0 

0.0500 

0.0554 

-0.0054 

0.9276 

372.0 

0.5300 

0.6705 

-0.1405 

0.8703 

372.0 

1.4200 

1.3938 

0.0262 

0.7974 

372.0 

3.2300 

3.0086 

0.2214 

0.7402 

372.0 

5.4700 

5.2162 

0.2538 

0.7194 

372.0 

6.5300 

6.3315 

0.1985 

0.7022 

372.0 

7.5900 

7.4222 

0.1678 

0.6859 

372.0 

8.8000 

8.6300 

0.1700 

0.6644 

372.0 

10.4800 

10.5047 

-0.0247 

FITTING  STATISTICS  FOR  5  PARAMETERS  FITTED  TO  55  DATA: 
Ru=  2.8227%  Rw=  4.5488% 

CHIA2  (NON-WEIGHTED)=  0.0347  CHIA2  (WEIGHTED)^  121.9185 
SUM(NON-W  DIFFSA2)=  1.73353  SUM(W-DIFFSA2)=  6095.9230 
MAXIMUM  DELTA-PRESSURE=-0.862 
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Table  Y.  The  Mie-Gruneisen  fitting  results:  Ko=A.9  GPa  (room-T),  Ko,=8.2  (room-T), 
Fo=1.005  (cm3/g),  Debye  temperature  6^0 =150  K  (fix),  q=l  (fix),  Gruneisen  parameter 
7=0.17  ±0.10.  In  table:  V  is  observed  volume,  P  is  observed  pressure,  Pcalc  is  calculated 
pressure,  AP=P-Pcalc,  and  APthermal=thermal  pressure. 

V(cm3/g)  T(K)  P  GPa  P(calc),  AP  APThermai(calc) 


0.992 

271. 

0.00 

-0.05 

0.05 

-0.07 

0.975 

271. 

0.05 

0.05 

0.00 

-0.12 

0.907 

271. 

0.53 

0.65 

-0.12 

-0.23 

0.855 

271. 

1.42 

1.38 

0.04 

-0.08 

0.784 

271. 

3.23 

3.11 

0.12 

0.00 

0.730 

271. 

5.47 

5.34 

0.13 

0.01 

0.711 

271. 

6.53 

6.41 

0.12 

0.00 

0.694 

271. 

7.59 

7.47 

0.12 

0.01 

0.677 

271. 

8.80 

8.73 

0.07 

-0.05 

0.657 

271. 

10.48 

10.55 

-0.07 

-0.19 

1.005 

302. 

0.00 

0.01 

-0.01 

0.00 

0.985 

302. 

0.05 

0.11 

-0.06 

-0.06 

0.911 

302. 

0.53 

0.73 

-0.20 

-0.19 

0.858 

302. 

1.42 

1.46 

-0.04 

-0.03 

0.787 

302. 

3.23 

3.14 

0.09 

0.10 

0.732 

302. 

5.47 

5.35 

0.12 

0.12 

0.712 

302. 

6.53 

6.44 

0.09 

0.10 

0.695 

302. 

7.59 

7.52 

0.07 

0.08 

0.679 

302. 

8.80 

8.75 

0.05 

0.06 

0.659 

302. 

10.48 

10.50 

-0.02 

-0.02 

1.011 

320. 

0.00 

0.05 

-0.05 

0.03 

0.990 

320. 

0.05 

0.16 

-0.11 

-0.03 

0.912 

320. 

0.53 

0.78 

-0.25 

-0.17 

0.858 

320. 

1.42 

1.53 

-0.11 

-0.03 

0.788 

320. 

3.23 

3.19 

0.04 

0.12 

0.736 

320. 

5.47 

5.24 

0.23 

0.31 

0.716 

320. 

6.53 

6.28 

0.25 

0.33 

0.699 

320. 

7.59 

7.33 

0.26 

0.34 

0.683 

320. 

8.80 

8.53 

0.27 

0.35 

0.662 

320. 

10.48 

10.25 

0.23 

0.31 

1.015 

332. 

0.00 

0.08 

-0.08 

0.05 

0.993 

332. 

0.05 

0.19 

-0.14 

-0.01 

0.912 

332. 

0.53 

0.83 

-0.30 

-0.17 

0.858 

332. 

1.42 

1.59 

-0.17 

-0.04 

0.787 

332. 

3.23 

3.28 

-0.05 

0.08 

0.733 

332. 

5.47 

5.45 

0.02 

0.15 

0.713 

332. 

6.53 

6.53 

0.00 

0.13 

0.695 

332. 

7.59 

7.64 

-0.05 

0.08 

0.679 

332. 

8.80 

8.83 

-0.03 

0.10 

0.655 

332. 

10.48 

10.96 

-0.48 

-0.35 

1.046 

372. 

0.00 

0.13 

-0.13 

0.17 

1.019 

372. 

0.05 

0.23 

-0.18 

0.12 

0.928 

372. 

0.53 

0.84 

-0.31 

-0.01 
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0.870 

372. 

1.42 

1.55 

-0.13 

0.17 

0.797 

372. 

3.23 

3.12 

0.11 

0.40 

0.740 

372. 

5.47 

5.25 

0.22 

0.52 

0.719 

372. 

6.53 

6.31 

0.22 

0.51 

0.702 

372. 

7.59 

7.35 

0.24 

0.53 

0.686 

372. 

8.80 

8.49 

0.31 

0.60 

0.664 

372. 

10.48 

10.27 

0.21 

0.51 

STANDARD  DEVIATION:  0.162514 
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